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j— H , We demonstrate that the space-time statistics of the birth of turbulent spots in boundary layers can 
' pe reconstructed qualitatively from the average behavior of macroscopic measures in the transition 
(Xione. The conclusion in [1] that there exists a connection between the patterns in laminar instability 
O^nd the birth of turbulent spots is strengthened. We examine why the relationship between instability 

* ^^nd transition to turbulence is manifest in some cases and appears to be totally absent in others. Novel 
^cellular automaton type simulations of the transition zone are conducted, and the pattern of spot 
j^Dirth is obtained from secondary instability analysis. The validity of the hypothesis of concentrated 

'"^breakdown, according to which most turbulent spots originate at a particular streamwise location, 
assessed. The predictions made lend themselves to straightforward experimental verification. 

<n : 
> ■ 

Introduction 

<n : 

^}rhe physics of the transition to turbulence in a boundary layer has been the 
(^subject of a large number of theoretical and experimental investigations. While 
l Omost of the process is understood, as shown in the schematic picture in fig- 
^^ure n there are still gaps in our knowledge, as we shall discuss below. The 

^process begins with the linear amplification of two-dimensional disturbance 
•^wfaves of a narrow band of frequencies. Once the disturbance modes have 

£>grown to a certain amplitude, they destabilize three-dimensional secondary 
F^modes [2-4]. These modes display a locally regular pattern of maxima, which 

ftftay be aligned (harmonic) or staggered (subharmonic) . Further downstream, 

£*im what is termed the transition zone, turbulent spots, or concentrated patches 
^o£ turbulence, are found, surrounded by quiet laminar flow. The spots grow as 

5-nthe convect downstream, and merge with each other to make the flow asymp- 
. .^totically fully turbulent. Our interest here is in the breakdown of instability 
waves into turbulent spots. We do not explain the complex physics occurring 
between the zone of secondary instability and the region containing spots, but 
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Figure 1. Sequence of events in the laminar-turbulent transition process, on a boundary layer 
formed by the flow past a semi-infinite plate. 



we present evidence showing that the pattern existing in the first is retained 
for the most part in the second. 

It is at present almost impossible either experimentally or numerically to 
track the birth and downstream growth of a large number of individual tur- 
bulent spots and to obtain detailed statistics. It is however not difficult to 
measure average quantities such as turbulent intermittency (the fraction of 
time that a flow is turbulent). With these, and knowing the behaviour of a 
typical spot, we demonstrate that it is possible to reconstruct the scenario of 
spot birth. It was recently shown [1] that intermittency behaviour in many 
flows is consistent with spot birth as dictated by the most unstable secondary 
mode. We present another example of this in decelerating flows. Hitherto, it 
was assumed that spot birth is completely random. No contradiction to this as- 
sumption was presented by intermittency measurements in low external noise 
environments, in the flow past a flat plate aligned with the flow. Since most 
intermittency measurements were made for this specific flow, the assumption 
of random spot birth gained credence and was not re-examined. We show that 
in this special case, the intermittency variation obtained from both random 
and regular spot birth are similar, and discuss why. We then present other 
measures which even in this case would differentiate clearly between random 
and regular spot birth. We exploit the versatility of our simulation method to 
demonstrate (section [7J) that the effect of geometry (in the form of transverse 
curvature) is significant in the boundary layers on cylinders. It is hoped that 
the present predictions will motivate experimental efforts. 

We conduct secondary instability computations in the standard manner (sec- 
tion EJ), and employ novel cellular automaton type stochastic simulations for 
spot birth, growth and convection, as described in section G3 We make a num- 
ber of assumptions. First, while the secondary instability analysis is three- 
dimensional, the stochastic simulations are conducted in a plane parallel to 
the solid body, when viewed from above. This means that the variations within 
the boundary layer in a direction y normal to the wall are averaged over. The 
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x-y profile of turbulent spots is not uniform, but given that they are flat objects 
and the boundary layer thickness is very small (0(R~ 1 ) compared to the other 
dimensions), averaging over y is a valid approximation. Secondly, we employ 
the hypothesis of concentrated breakdown, according to which all spots are 
born at a particular streamwise location [5], i.e. at a particular Reynolds num- 
ber. We show in section El that this is a reasonable assumption. We also make 
use of the experimentally and numerically observed shape and downstream 
behaviour of single spots, especially the fact that spot growth is self-similar at 
any pressure gradient. 

A spatial pattern in turbulent in the form of spirals and spots are observed in 
the experiments of Prigent et al. [6] in plane Couette flow and Taylor-Couette 
flow. They found that the spatial modulation of turbulent intensity obeys 
the dynamics of coupled amplitude (Ginzburg-Landau type) equations with 
noise. It is important to note that the simulations of [7, 8] in boundary layers 
establish that structures in the unstable laminar region are the precursors of 
turbulent spots. These simulations were of bypass transition. While the present 
work is about the traditional Tollmien-Schlichting mechanism taking place in 
low-disturbance environments, the basic connection between instability and 
transition is similar. 



2 Analysis of secondary instability 

The growth of small (linear) perturbations in a boundary layer, under the 
assumption of parallel flow, is described by the Orr-Sommerfeld equation [9,10] 



where U(y) is the mean velocity and R the Reynolds number, = U !X (x)9(x)/u, 
Uqc is the local freestream velocity, 6 is a characteristic boundary layer thick- 
ness, and v the kinematic viscosity. v(y) is the amplitude of the normal com- 
ponent v of the disturbance velocity, given by 



k and oj being the wavenumber and the frequency of the disturbance mode 
being considered. The primes denote differentiation with respect to y. 

When the linear modes have grown to a significant amplitude (of the order 
of a percent of the mean flow), the new flow becomes unstable to secondary, 
typically three-dimensional, disturbances. The linearised equations describing 
the secondary instability are obtained by writing flow quantities in the form 




(1) 



v(x, y, t) = v(y) exp[i(kx — u>t)] 
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[11,12] 



U(x,y,t) = { U(y)z + A u(y)i + v{y)j exp [1 (kx - ut)] \ +u s (x, y, z, t). (2) 



The quantity in the curly brackets is the basic flow consisting of the laminar 
profile plus the primary instability wave. The primary instability is very slowly 
growing, and the flow may be taken to be periodic. A stands for the amplitude 
of the primary disturbance. The three-dimensional secondary disturbance is 
written as 



u s (y, r ± , t) = Re\ u s+ (y) exp [1 (k+ • r± - u + t)]+u s -(y) exp [1 (k_ • r± - u-t)] 



into the Navier-Stokes and continuity equations, retain linear terms in the 
secondary, and eliminate the disturbance pressure and streamwise component 
of the velocity. On averaging over x, z and t, only the resonant modes survive, 
which are related as follows: 



where rj_ 




k + + k- = k, 



and 



(4) 



We then obtain 



1 



(uj + - k + U) + ]-{D 2 -k\- [(k 2 + + p 2 ) f s+ - Dv s+ ] - ik+U'v s+ 




[\k + uD + vD 2 + ik-Du] v*_ + (/3 2 - k^k + ) vD+ 





and ^ - k+ U) + ±-(D 2 - k 2 + - p 2 ) (Df s+ -v s+ )-ik + U'f s+ 
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The quantity f s +(= —^w s+ /f3) is proportional to the spanwise component 
of the secondary disturbance velocity, and the operator D stands for differen- 
tiation with respect to y. The boundary conditions are 



Us 



at y = & u s — > as y — > oo. 



(7) 



Equations © and (jHJ), along with two corresponding equations in v*_ and 
/*, describe an eigenvalue problem for the secondary instability. Note that 
secondary instabilities account for one contribution to nonlinearity, that of 
the triad interaction between the primary Tollmien-Schlichting mode and two 
three-dimensional disturbances, the sum of whose wave-numbers equals that 
of the primary. 

A Chebyshev spectral collocation method is used for the numerical solution 
[13]. A grid stretching given by 



Vj 



b(l + Y L 



1 + 2b/ Voo - Yj 



(8) 



is applied, where 



Yj = cos(-^-) 



N„ 



0,1, 2,..., N s 



9-' 



(9) 



are Chebyshev collocation points. We thus transform the computational do- 
main from (—1,1) to (0, yoo), where yoo is chosen to be at least 5 times the 
boundary layer thickness, and cluster grid points close to the wall by tuning 
the parameter b. Results for the boundary layer over a flat-plate compare very 
well with those of [11]. On changing the number of grid points N g from 80 to 
160, the eigenvalues remained identical up to the sixth decimal place. 

We are interested here in the flow past flat solid surfaces inclined at an angle 
2m/ (m + 1) to the flow, where the external flow velocity varies downstream 
as C/qo ~ x m [14]. We consider constant pressure (m = 0) and decelerating 
boundary layers (m < 0) in this paper, since the behaviour in the transition 
zone differs in the two. The profile of mean streamwise velocity for a given m 
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is described by the Falkner-Skan equation [15] 

r+/r + ^ T (i-/ /2 ) = o (io) 

where U = /'. The boundary conditions are / = /' = at the wall, and 
/' — > 1 as y — > oo. Computations of instability have been carried out for 
a variety of Reynolds number and primary disturbance amplitude. The sub- 
harmonic mode, of wavenumber (= k + = k-) = k/2 is found always to be the 
most unstable. 



3 The transition zone 

The transition zone is easiest to describe quantitatively in terms of the vari- 
ation with the streamwise coordinate, x of the intermittency, 7. Given that 
the Reynolds number R(x) is the only parameter in the instability problem, 
it is reasonable to assume that most turbulent spots will originate within a 
narrow spanwise strip around a particular streamwise location xt, denoted 
here as the location of transition onset. This is the hypothesis of concentrated 
breakdown [5]. With spots appearing randomly in accordance with a Poisson 
distribution in time and a uniform distribution in the spanwise coordinate z, 
7 would vary downstream as [5] 



7 = 1 — exp 



na 2 
(x - x t ) 



(11) 



or 



F = v /-log(l- 7 ) = 



— na 



(x - x t ) 



(12) 



such that the intermittency parameter F varies linearly in x. Here Uoo is the 
free stream velocity, n is the number of spots forming per unit length in the 
spanwise direction, z, per unit time, at x t . The shape and downstream growth 
of turbulent spots is experimentally observed to be similar in all pressure 
gradients: a turbulent spot maintains an arrowhead shape when viewed from 
above, and remains self-similar as it grows [16,17]. It is convenient therefore 
to define a non-dimensional spot propagation parameter as 



a = 



' 1 

U~ r 



1 " 



Uoo tan(, 



(13) 
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where ( is the angle subtended by the spot at its origin and Uh and U r respec- 
tively are the speeds with which its head and rear convect. 

It is relatively straightforward to measure other average transition zone 
quantities such as the burst rate and the persistence time distribution of lam- 
inar or turbulent flow as functions of the streamwise distance. The burst rate 
B is here defined as the mean number of switches from laminar to turbulent 
flow per unit time at a given location, and the persistence time W of laminar 
flow is the duration of time for which the flow remains continuously laminar at 
a given location. The transitional intermittency has been measured in many 
experiments, but there are no data available for the other measures, as far as 
we know. It is our suggestion that a lot of information can be gained from 
knowing how these quantities vary in the streamwise direction. For a Poisson 
process, the burst rate is related to the local intermittency by 



and the probability density function of persistence time W of laminar flow is 



In quiet facilities, the parameter F in equation El in constant pressure flow 
(to be discussed later) varies almost perfectly linearly with x, but its variation 
in decelerating flows is nonlinear with an increasing slope, as seen in figure El 
The explanation [1] for the nonlinear behaviour of intermittency is that spot 
birth is not random, but mainly regular, in a pattern dictated by the most 
unstable secondary mode. 



4 Stochastic Simulations 

Stochastic simulations of the transition zone (using a cellular automaton-like 
approach) are performed, employing observed properties of spot growth and 
convection [1,19]. We consider two scenarios for spot birth. The first is a ran- 
dom breakdown, as discussed above, with no relation to instability. The second 
is predominantly regular breakdown, where the pattern prescribed is that of 
the maxima in disturbance vorticity in the secondary instability. Simulated 
transition zones for random and regular (harmonic) breakdown are shown 
schematically figures 01 and |U respectively. 

The zone, as viewed from above, is discretized into 200 x 400 rectangular 
grids in x and z. Each grid is assigned an integer x, equal to if the flow 
there is laminar and 1 if it is turbulent. A regular breakdown scenario is 
prescribed as follows. During the first time interval, one spot is generated 



£<x(l- 7 )(-ln(l- 7 )) 



1/2 



(14) 




(15) 



February 2, 2008 5:31 Journal of Turbulence vinod 



2.5 p 

2- 

1.5 - 
F 

1 

0.5 



0.2 



0.4 0.6 0.8 



Figure 2. Variation of the intermittency parameter F in adverse pressure gradient, m - 

Symbols: experiment [18]. 



-0.05. 



every N z sites, at X = 1 and Z = IN Z + i, I = 0, 1, 2 • • ■ L z /N z , where i is an 
integer between 1 and N z . The uppercase stands for discretised coordinates. 
After N t time steps the spots form at spanwise locations which are staggered 
in the spanwise coordinate with respect to the spots generated earlier, i.e., at 
X = 1 and Z = lN z +j for harmonic and Z = (1 + 1/2)N Z +j for subharmonic 
breakdown where j = i + 1. This prescription corresponds to spots forming at 
the crests of an oblique wave of spanwise and streamwise wavenumbers 



2vr 



N Z AZ 



and 



k = 



2vr 



N t ATv 



(16) 



respectively, where v is the streamwise velocity of the wave crest. The phase 
of the oblique wave is randomised with a small probability (1%), and a small 
fraction (5%) of the spots are generated randomly. It is assumed that this pre- 
scription mimics qualitatively the randomisation due to external disturbance. 
Changes in these fractions do not change the answers significantly. 

In a random breakdown, spots appear at X = 1 in accordance with a Pois- 
son distribution in time and a uniform distribution in the spanwise direction. 
In each case, the downstream growth of turbulent spots follows the same al- 
gorithm. During one time interval, the front of each spot moves ahead by 2 
grid locations, while the rear moves forward by 1 grid location. The lateral 
dimension of the spot is increased by 1 grid location on either side. The simu- 
lated spot is thus triangular, and retains its shape for all time. A real spot is a 
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Figure 3. Schematic diagram of the simulation domain in a random spot breakdown scenario. 
Spots appear according to a Poisson distribution in time, and are uniformly distributed in the 

spanwise direction. 



Figure 4. Schematic diagram of the simulation domain in regular breakdown. Spots are formed at 
regular intervals in the spanwise direction and are staggered after every time period, in a manner 

dictated by the secondary instability. 



blunter at its leading edge than a triangle, but this is assumed not to change 
the results significantly. 

In experiments, the speed of the front of the spot is comparable to the flow 
velocity t/oo while the rear moves forward approximately at 0.5 [Too. The half- 
angle ( subtended by the spot at its origin is about 10° in a constant pressure 
boundary layer, and increases with adverse pressure gradient. 

In Figs. El and El it is demonstrated that intermittency behaviour is sensitive 
to the ratio of spanwise to streamwise wavenumbers. This observation is used 
to strengthen the connection we make between instability and its signature in 
the transition zone. The simulations shown in Fig. display a distinct change 
in slope at X = 13, just downstream of the location {X = 11) at which the 
heads of a set of spots formed at a time T just touch the rears of the spots 
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x 

Figure 5. Streamwise variation of the intermittency parameter for N z = 49 and Nt = 10. Symbols: 
regular breakdown, solid line: random breakdown. The mean spot generation rate is kept constant 
at N = 0.01. The dashed line shows the best linear fit for the downstream part of F. 



formed at T — NfAT. In the simulations of Fig. H3 all parameters have been 
kept the same, except that the values of N z and Nt have been switched. In 
this case the spots touch each other laterally first (at X = 6) downstream 
of which a sharp change in slope is again evident. At low F, i.e., upstream 
of spot merger, the intermittency behaviour in each case is exactly the same 
as that when spots form randomly (with the same breakdown rate). This is 
not surprising, given that the spots are too small to "see" each other. At 
higher intermittencies, quite surprisingly, for the same mean breakdown rate 
and identical spot growth, the length of the transition zone is seen to be very 
different, which tells us that spot-merger plays a large role in transition. A 
purely lateral merger results in a much higher degree of overlap (on both sides), 
and consequently, less of the region is occupied by spots. This of course means 
that transition to turbulence will proceed much more slowly as shown in Fig. [HJ 
A purely longitudinal merger on the other hand, results in far less spot overlap, 
meaning more occupancy by spots, and a shortened transition zone (Fig. 
Incidentally, it is often the case that experimental results seem to lie on two 
straight lines, as in this figure. This has been termed a 'sub-transition' [20]. 
Our simulations show that this change in slope in not a result of spots suddenly 
changing their downstream growth behaviour, but is a natural consequence of 
their spatial arrangement. 

We now ask which kind of merger scenario would be consistent with growing 
instability modes. To do this, we perform a secondary instability analysis to ob- 
tain the streamwise and spanwise wavenumbers for which the secondary wave 
grows fastest. For example, for an adverse pressure gradient flow of Falkner- 
Skan parameter m = —0.06, the maximum growth rate of the secondary 
instability wave is 0.067, and occurs when the wavenumbers (k and (3) are 
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Figure 6. Streamwise variation of the intermittency parameter for N z = 10 and Nt = 49. Other 
parameters and notation are the same as in figure [S] It can be seen that the transition zone is much 

longer in this case. 



0.165 and 0.215 respectively. The corresponding phase speed is 0.524[/oo. Us- 
ing the fact that the rear of a spot travels at approximately half the freestream 
velocity, we have tan£ = 2AZ /U^AT '. The spot propagation angle is taken as 
~ 20°, which is close to the experimental value for this pressure gradient [21], 
and from equation ()16|) we may estimate the breakdown ratio N z /N t to be 
~ 4.9. This ratio changes only by 2% for a 50% change in Reynolds number, 
and may be taken to be independent of Reynolds number. 

In figure [7| we consider an adverse pressure gradient, m = —0.05. The ra- 
tio N z /Nt obtained from stability analysis for this case is 3.23, and in the 
stochastic simulations we use N z = 39 and Nt = 12. The analysis is done at a 
Reynolds number of 220, based on the boundary layer momentum thickness, 
to match the experimental value. A similar result for m = —0.06 is available 
in [1]. The results of both simulations are in good agreement with experiments. 
In particular, the changing slope of the intermittency parameter F(x) is fol- 
lowed well. We may conclude that a spot breakdown pattern as prescribed by 
the secondary instability gives rise to intermittency behaviour in qualitative 
agreement with measurements in highly decelerating flows. 

It may be noticed that there is some discrepancy at the beginning of the 
transition zone, where it is known that experimental data may be inaccurate: 
errors could arise unless the data is collected over extremely long times. Sec- 
ondly, tiny patches of turbulence can be difficult to distinguish from noise. It 
is also to be remembered that the hypothesis of concentrated breakdown is 
an idealisation: in reality spots would be forming within a narrow streamwise 
strip around x = Xt, rather than at one particular location, this can smear out 
the intermittency at 7 ~ 0. 

The burst rates obtained from the present simulations with a random break- 
down are in agreement with equation 1141 as seen in figure |H1 Any deviation 
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from this behaviour is a sign that spot breakdown is not random. When the 
breakdown is mostly regular, it is clear that the burst rate distribution is 
"peaky" and symmetric with respect to 7 = 0.5. The burst rate variation with 
intermittency is similar for sub-harmonic and harmonic type breakdowns. We 
suggest here that there is a need for experimental measures of these quanti- 
ties, especially in adverse pressure gradient boundary layers, which will help 
us quantify the breakdown scenario. The burst rate behaviour displayed when 
the dominant merger is lateral is qualitatively different from the case where the 
merger is longitudinal. The burst rates obtained from the simulations whose 
intermittencies are shown in figure |S] and H3 are plotted in figure El and are a 
demonstration of this. 

The probability density function of the persistence time W is plotted in 
figure IH1 The data is obtained by running the simulation over 20 million time 
steps (after reaching stationary state), monitoring a particular streamwise lo- 
cation in the middle of the span and collecting statistics of lengths of strings of 
zeroes between two l's. The reverse information, about strings of l's between 
two zeroes, is less interesting. We present results for W at location where the 
intermittency is 7 = 0.1. For a random spot breakdown, the probability den- 
sity of the persistence time decays exponentially with persistence time, as is 
expected from equation (|15|) . At low intermittency levels, in the case of regular 
breakdown, an overall decay is significantly modulated by ups and downs. At 
higher levels of intermittency, if the breakdown is predominantly regular, the 
probability of very large waiting times is extremely low, as expected. The most 
probable persistence times correspond to the modulated streamwise extent of 
the laminar zone between two rows of spots. In the case of subharmonic (stag- 
gered) spot arrangement, there is a probability of skipping one row of spots, 
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Figure 8. Variation of the burst rate with intermittency. The solid line shows a subharmonic 
breakdown, with Nz/Nt = 4.9, with 5% randomness. The scale on the ordinate is arbitrary. The 
result is similar to the predominantly regular harmonic breakdown case shown in [1]. The dashed 
line shows equation 1141 and the circles are simulation results with random breakdown. 



10 2 F 




Figure 9. Probability density function of persistence time of laminar flow at (a) 7 = 0.1. Dashed 
line: random breakdown, solid line: 90% regular (sub-harmonic), long dashes: 50% regular. 

giving larger persistence times. In figure E3 the probability density function 
of the persistence time from a simulation where the lateral merger is dom- 
inant is plotted. The qualitative behaviour is different and an experimental 
measurement of p(W) is can thus indicate the dominant breakdown pattern. 
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Figure 10. Probability density function of persistence time of laminar flow at 7 = 0.1 (solid line) 
and 7 = 0.5 (dashed line), The simulations are carried out under conditions where lateral merger is 

dominant as in Fig. l6l 



5 Is spot birth random in constant pressure gradient flows? 

In constant pressure (flat plate boundary layer) flow, two kinds of behaviour 
are observed. In high-disturbance environments, the intermittency behaviour 
is similar to that seen in adverse pressure gradients, except the dominant 
mechanism is the secondary instability of streaks. The experimentally mea- 
sured intermittency in such a flow [23] has been shown [1] to be consistent 
with the dominant instability pattern. Most experiments, however, are con- 
ducted in flat plate in environments which are painstakingly maintained quiet. 
Here, the intermittency parameter F is linear in x, consistent with random 
spot-birth. Do secondary instabilities not directly trigger spot birth here? A 
partial answer may be gleaned from figure ITTl where the spanwise/ stream- 
wise ratio of wavelengths of the dominant mode is shown as a function of 
the pressure gradient. The ratio decreases with decreasing pressure gradient, 
and N z /Nt is about 2.3 for a constant pressure boundary layer (since now 
k = 0.085, (3 = 0.146 and v = 0.353 for the most dangerous mode). For a 
regular breakdown with this ratio, simulations show that the variation in F is 
linear up to about 7 = 0.75 ( figure IT2*|) . the slope increases towards the end 
of the transition zone. On the lines of the discussion on Fig. Inland |SJ we may 
expect that spot mergers are now a combination of lateral and longitudinal, 
with longitudinal being marginally dominant. The weakened dominance means 
that the slope change occurs much closer to full turbulence, at 7 ~ 0.8. 

The obliqueness of the most unstable wave provides only a partial answer 
for why the variation of F is linear in a constant pressure boundary layer. 
The perfectly linear fit in a vast number of experimental data is also likely to 
be because the degree of randomness in this case is actually higher than in 
an adverse pressure gradient boundary layer, as suggested by Gostelow and 




February 2, 2008 5:31 Journal of Turbulence vinod 



15 



N 

Z 

N. 



-0.06 -0.05 -0.04 -0.03 -0.02 -0.01 

m 

Figure 11. The variation with pressure gradient of the ratio N z /Nt as determined by the most 

unstable secondary mode. 
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Figure 12. Variation of intermittency parameter F in the boundary layer over a flat plate. Solid 
line: stochastic simulations with N z = 32, Nt = 13 as dictated by secondary instability, the dots are 
experimental data in constant pressure boundary layers [24]. 



Walker [25] . The following qualitative arguments indicate why we may expect 
this. 

A decelerating flow is inherently much more unstable [26, 27], as demon- 
strated in figures an d 1141 Since the primary wave is slowly growing, we 
assume that the secondary mode at a given time can be computed by taking 
the primary wave to be of constant amplitude (this is a temporal analog of 
the "parallel- flow" assumption often employed in spatially-developing flow). 
It is to be noted that the growth of the secondary mode can be faster than 
exponential if the primary mode is unstable as well. Within a very short down- 
stream distance, the disturbance waves in an adverse pressure gradient flow 
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Figure 13. Growth rate of most dangerous secondary sub-harmonic mode at a Reynolds number 
R = 600, based on the momentum thickness of the boundary layer and Uoo ■ Solid line: 
m = -0.06, k+ = 0.14, dashed line: m = 0, fc+ = 0.085. 



Figure 14. The variation of amplitude of secondary disturbance wave with time in adverse 
(m = -0.06, k = 0.185 and (3 = 0.12,solid line) and zero (k = 0.13 and f3 = 0.805, dashed line) 

pressure gradient boundary layers. 

achieve the threshold amplitude required for breaking down into spots. In the 
constant pressure case the attainment of the required threshold is much slower, 
offering greater opportunity for stochastic effects. 



A significant amplitude modulation of wave packets could be another cause 




350 



time 



February 2, 2008 5:31 Journal of Turbulence vinod 



17 



Pressure gradient (m) 


Reynolds number 


phase speed 


group velocity 
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Tabic 1. Typical phase and group velocities in zero pressure-gradient and decelerating boundary layers 
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Tabic 2. Effect of the primary disturbance amplitude on the obliqueness of the most dangerous mode. 

for randomness in zero pressure gradient flow. Amplitude modulation would 
result in different waves reaching the threshold at different streamwise stations, 
and spot birth would be smeared out over a streamwise distance. In the case 
of decelerating flow, the group velocities of the waves are much larger, as seen 
in Table I, and the amplitude modulation takes place over a much narrower 
width. At a Reynolds number of 200, the group velocity c g of the most unstable 
secondary in adverse pressure gradient is more than twice the corresponding 
wave over a flat plate. At a Reynolds number of 600, this ratio is more than 
3. The phase speeds on the other hand are not that different. 

Finally, the wavelength ratio of the dominant mode varies significantly with 
primary disturbance amplitude in flat-plate flow (see [28]). This, however, is 
not the case in strongly decelerating flow, as seen in Table II. If the external 
noise is irregular in amplitude, a greater randomisation would take place in 
flat plate flow than in decelerating flow. 

These results indicate that in a decelerating boundary layer, formed by the 
flow past an inclined plate, the connection between instability and transition is 
likely to be much easier to observe, and we therefore recommend experimental 
work in adverse pressure gradient, rather than on flat plate boundary layers. 

6 Effect of concentrated breakdown 

Narasimha [5] proposed the hypothesis of concentrated breakdown in which 
all spots form within a narrow spanwise belt around the location of transition 
onset, i.e., within xt ± e, where e << xt- While the intermittency behaviour 
resulting from the hypothesis and the assumption of a random breakdown 
matches experiment very well in constant pressure and low disturbance en- 
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vironments, its validity has been a matter for debate, see e.g. [29]. Since a 
particular disturbance amplitude is achieved at a particular Reynolds num- 
ber, it is plausible at least that upstream of a given x location, no spots will 
form. A large number of spots are likely to form in the vicinity of xt, but 
given that instability modes which are not dominant could continue to grow 
in the transition zone, a small number could be born at any location down- 
stream. We conduct simulations to estimate how much the intermittency and 
other parameters depend on making this hypothesis. We do this by allowing 
an increasing fraction of the spots to form with equal probability anywhere 
downstream of xt- A random breakdown is prescribed. Our conclusion is that 
even in the unlikely event of a large fraction of spots being born downstream, 
the intermittency is dominated by spots forming at xt, since at any given x, 
these spots are much larger than those born more recently. The intermittency 
distributions for the extreme cases tried are shown in figure ITol It is seen that 
even when 80% spots form downstream of the onset location, the variation of 
F is still practically linear. However when all the spots are born downstream 
of xt (as modelled by Emmons [30]), there is a noticeable departure from lin- 
earity. In an instability-driven transition, it is highly likely, as discussed above, 
that a significant fraction of spots will be born around xt- It may be noticed 
that the intermittency behaviour in the case of Emmons breakdown does not 
appear very different qualitatively from a regular breakdown, and would thus 
be difficult to distinguish in an experiment. However, an experiment which 
measures persistence times can distinguish very simply between them, as is 
evident from figure [TBI In an Emmons breakdown, the persistence time distri- 
bution decays exponentially, but with a smaller slope than for the concentrated 
breakdown, as expected. 



7 Axisymmetric Boundary Layers 

Transition to turbulence in the boundary layers forming around axisymmetric 
bodies is poorly understood in spite of wide application in the motion of 
submarines, fishes etc. When the transverse curvature is significant, transition 
can proceed quite differently from a two-dimensional boundary layer [20,31]. 
When the typical patch of turbulence, consisting either of a single spot or a 
group of spots which have merged laterally, attains a width of the order of the 
diameter of the cylinder, it wraps itself around the body. Downstream of the 
location of wrap, further lateral growth is not possible. The turbulent patch 
then resembles a sleeve [32] displaying only a one-dimensional growth in the 
streamwise direction. 

We have carried out stochastic simulations of the birth and downstream 
propagation and growth of turbulent spots in the transition zone of an ax- 
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Figure 15. Effect of the hypothesis of concentrated breakdown on the intermittency distribution. 
Solid line: spots are allowed to form only at xt- Dashed line: 80% spots are born downstream of the 
onset location. Dot-dashed line: all spots form downstream of the transition onset (Emmons 

breakdown) . 




W 

Figure 16. Effect of concentrated breakdown on persistence time distribution. The persistence 
time is computed at location where intermittency 7 = 0.1. The curve marked R2 is according to the 

hypothesis of concentrated breakdown, Rl is for random breakdown anywhere downstream of 
xt [30], PI represents periodic breakdown (sub- harmonic), the pattern is obtained from secondary 

instability. 

isymmetric boundary layer. The downstream variation of the intermittency 
parameter F is shown in figure 1171 The quantity c is the circumference of 
the body. In the initial region, transition proceeds exactly as it would in two- 
dimensional flow. This is because spots are too small to "see" the body. When 
spots wrap themselves around the cylinder, a qualitative change is observed, 
as is expected from the discussion above, and transition proceeds much more 
slowly after this. The burst rate B shown in figure^] is another indication of 
the differences in the transition process. 
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Figure 17. The intermittency factor F Vs. x, for different circumferences of the cylinder. The 
straight line is the result of two-dimensional simulations with the same spot birth rate. 

8 Summary and discussion 

Stochastic simulations, inspired by a cellular-automaton approach, of the gen- 
eration and propagation of turbulent spots in transitional boundary layers have 
been conducted, employing the observation that spot growth is self-similar. It 
is demonstrated that the pattern of spot birth may be inferred from the down- 
stream variation of average quantities in the transition zone, such as the inter- 
mittency and the burst rate. This is because the qualitative behaviour in the 
transition zone depends on whether there is a pattern in spot merger, and if so, 
whether lateral or longitudinal merger is dominant. Contrary to present belief, 
our results indicate that relatively simple experiments can tell us a great deal 
about the connection between instability and transition to turbulence. We 
show that experiments conducted in decelerating flows are much more con- 
ducive to exploring this connection than constant pressure flows. Transverse 
curvature has the effect of slowing down the transition process. 

The simulations of [7] showed a degree of spanwise periodicity in the ar- 
rangement of the backward jets (the precursors to spots). Only some of the 
jets however gave rise to turbulent spots. It has been our contention in TS- 
mode driven transition that the degree of randomness in the birth of spots is 
expected to decrease as the pressure gradient becomes increasingly adverse. It 
would be interesting to see whether this is true of bypass transition as well. A 
direct numerical simulation of an adverse pressure gradient flow should, if this 
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Figure 18. The burst rate in the transition zone of the boundary layer around a cylinder. 

contention is right, produce turbulent spots from many more of the reverse 
jets. 
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